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Abstract 

The generation of hydrodynamic radiation in interactions of pulsed proton and laser beams with matter is explored. The beams 
were directed into a water target and the resulting acoustic signals were recorded with pressure sensitive sensors. Measurements 
were performed with varying pulse energies, sensor positions, beam diameters and temperatures. The obtained data are matched 
by simulation results based on the thermo-acoustic model with uncertainties at a level of 10%. The results imply that the primary 
mechanism for sound generation by the energy deposition of particles propagating in water is the local heating of the medium. The 
heating results in a fast expansion or contraction and a pressure pulse of bipolar shape is emitted into the surrounding medium. 
An interesting, widely discussed application of this effect could be the detection of ultra-high energetic cosmic neutrinos in future 
large-scale acoustic neutrino detectors. For this application a validation of the sound generation mechanism to high accuracy, as 
achieved with the experiments discussed in this article, is of high importance. 
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PACS: 


1. Introduction 

In 1957 G.A. Askaryan pointed out that ionisation and cavi¬ 
tation along a track of an ionising particle through a liquid leads 
to hydrodynamic radiation HI. In the 1960s, 1970s and 1980s, 
theoretical and experimental studies have been performed on 
the hydrodynamic radiation of beams and particles traversing 
dense media SH. 

The interest in characterising the properties of the acoustic 
radiation was, among other reasons, lead by the idea that the 
effect can be utilised to detect ultra-high energy (E > 10** eV) 
cosmic, i.e. astrophysical neutrinos, in dense media like wa¬ 
ter, ice and salt. In the 1970s this idea was discussed within 
the DUMAND optical neutrino detector project 0 and has 
been studied in connection with Cherenkov neutrino detector 
projects since. The detection of such neutrinos is consider¬ 
ably more challenging than the search for high-energy neutrinos 
{E > 10*** eV) as currently pursued by under-ice and under¬ 
water Cherenkov neutrino telescopes EMI- Due to the low 
expected fluxes, detector sizes exceeding 100 km* are needed 
Ha. However, the properties of the acoustic method allow for 
sparsely instrumented arrays with ~100 sensors/km*. 

To study the feasibility of a detection method based on 
acoustic signals it is necessary to understand the properties of 
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the sound generation by comparing measurements and simula¬ 
tions based on theoretical models. According to the so-called 
thermo-acoustic model 111, the energy deposition of particles 
traversing liquids leads to a local heating of the medium which 
can be regarded as instantaneous with respect to the hydrody¬ 
namic time scales. Due to the temperature change the medium 
expands or contracts according to its bulk volume expansion co¬ 
efficient a. The accelerated motion of the heated medium gen¬ 
erates an ultrasonic pulse whose temporal signature is bipolar 
and which propagates in the volume. Coherent superposition 
of the elementary sound waves, produced over the cylindrical 
volume of the energy deposition, leads to a propagation within 
a flat disk-like volume in the direction perpendicular to the axis 
of the particle shower. 

In this study, the hydrodynamic signal generation by two 
types of beams, interacting with a water target, was investi¬ 
gated; pulsed protons and a pulsed laser, mimicking the for¬ 
mation of a hadronic cascade from a neutrino interaction under 
laboratory conditions. With respect to the aforementioned ex¬ 
perimental studies of the thermo-acoustic model, the work pre¬ 
sented here can make use of previously unavailable advanced 
tools such as GEANT4 03 for the simulation of proton- 
induced hadronic showers in water. Good agreement was found 
in the comparison of the measured signal properties with the 
simulation results, providing confidence to apply similar sim¬ 
ulation methods in the context of acoustic detection of ultra- 
high energy neutrinos. A puzzling feature observed in previous 
studies—a non-vanishing signal amplitude at a temperature of 
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4°C, where for water at its highest density no thermo-acoustic 
signal should be present—was investigated in detail. Such a 
residual signal was also observed for the proton beam experi¬ 
ment described in this article, but not for the laser beam, indi¬ 
cating that the formation of this signal is related to the charge 
or the mass of the protons. 

2. Derivation of the Model 

In the following, the thermo-acoustic model EO is derived 
from basic assumptions, using a hydrodynamic approach. Basis 
is the momentum conservation, i.e., the Euler Equation 

Ot ^ OXj 

J=l J 

for mass density p, velocity vector field of the medium 
1 ^ 2 . V 3 ) and momentum-density tensor 

P^ij+pViVj ( 2 ) 


In the acoustic regime, where 6Q - 0, the momentum den¬ 
sity tensor can be expressed as 114 = pdij - v^pdij (us¬ 
ing Eqs. and where we assume an adiabatic density 
change with pressure. The non-linear kinetic term pvivj enter¬ 
ing nfj according to Eq. (j^ can be neglected for small devia¬ 
tions p' - p — pq from the static density po and small pressure 
differences p' - p - Po from the static pressure po 03. 

In the region B, where non-equilibrium deposition occur, one 
may make the ansatz 

Ylfj ^ p6ij + puiUj (7) 

with the direction m, of the beam which breaks the isotropy of 
the energy-momentum tensor and describes with the parameter 
P in an effective way the momentum transfer on the fluid. 

Although in non-equilibirum we apply Eq. with the en¬ 
ergy deposition density e = QjY of the beam. Then, with the 
additional energy-momentum tensor due to the beam 

( 511 ® := v] — e6ij + puiUj 


including the pressure p E). Equation Q can be derived from 
momentum conservation. In the derivation, energy dissipation 
resulting from processes such as internal friction or heat trans¬ 
fer are neglected. Motions described by the Euler Equation 
hence are adiabatic. Taking the three partial derivatives of Eq. 
Q with respect to Xj and using for the density the continuity 
equation 




d 


a non-linear wave equation can be derived: 

d^Uu 


^ = V 


df- 


. . , dxidxj 


arated spatial regions: Eirstly, a region B (‘beam’), where the 
energy is deposited in the beam interactions with the fluid and 
thus the wave excited in a non-equilibrium process; and sec¬ 
ondly, a hydrodynamic (‘acoustic’) region A, where the acous¬ 
tic wave propagates through the medium and where linear hy¬ 
drodynamics in local equilibrium can be assumed. This split¬ 
ting can be reflected by the momentum density tensor, rewriting 
it as 

for r € A 

(5) 


n,v(r^ - 




for r e B 

In local equilibrium the changes in mass density are given by 


the wave equation Q reads 

1 d^p' 
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The general solution for the wave equation can be written using 


a Green function approach as 

w ,, 1 r 

1 ^ik 

, 1 d^5W,(f,t') 

|r-r'| dx'.dx'i^ 

nm d^Snl(r^,n 

\r- r'l v^df^ 

(d6u^pf,t') dn® y 

\r- r'p 

^ Vsdt' ' |r-r'|j 


(9) 


with the components of the unit vector n, = (x, - xj)/|r - r'\ 
and the retarded time t' - t - \r- r'l/v^. Eor the last conver¬ 
sion, partial integration and the total derivative 
have been used repeatedly. Note that <511®. = 0 for r e A, so 
that the integration is carried out over the volume of the energy 
deposition region B. 

Assuming an energy deposition without momentum transfer 
to the medium, the kinetic term in the ansatz 0 can be ne¬ 
glected (J3 - 0) yielding 

p'(i^,t) = r 7^ ^ e (?',0 ( 10 ) 

Jy |r - r'l dP 


dp 


dp 


dp + 

S.N 


dS 



-^dp- 


a 6Q 


(6) 


with the bulk volume expansion coefficient a = ~p^\p,N^ the 
energy deposition 6Q - TdS, the adiabatic speed of sound Vj 
in the medium and the specihc heat Cp - ^%\p,N- 


for a thermo-acoustic wave generated solely by heating of the 
medium. The signal amplitude p' can be shown to be pro¬ 
portional to the dimensionless quantity v^ajcp when solving 
Eq. ( |T0) l for the case of an instantaneous energy deposition. 

Equation ( [T0| is equivalent to the results obtained from the 
approaches presented in ||2ll3- The derivation pursued above, 
however, uses a different approach starting with the Euler equa¬ 
tion and an anisotropic energy-momentum tensor, yielding a 
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more general expression in Eq. 0 . Only if assuming an 
isotropic energy deposition one arrives at the expression for the 
pressure deviation p' given in Eq. 

Note that the validation of the last assumption p - 0 being a 
good approximation would require a detailed knowledge of the 
momentum transfer from the beam to the medium. Taking it 
into account would result in an additional dipol term ~ 0 {u-n0 
in Eq. ( [T0| i which may become the dominant contribution to 
wave generation if a a; 0 close to 4.0°C. However, for 0-0 the 
pressure held resulting from a beam interaction in a medium is 
determined by the spatial and temporal distribution of the en¬ 
ergy deposition density e alone. The amplitude of the resulting 
acoustic wave is governed by the thermodynamic properties Vj, 
Cp and a, the latter three depending primarily on the tempera¬ 
ture of the medium. A controlled variation of these parameters 
in the conducted laboratory experiments and a study of the re¬ 
sulting pressure signals therefore allows for a precise test of the 
thermo-acoustic model. 

Simulations based on the thermo-acoustic model, as per¬ 
formed to interpret the results of the experiments described in 
the next section, will be discussed in Sec. Note that the en¬ 
ergy deposition density e and its temporal evolution for the pro¬ 
ton and laser beam interactions discussed in this paper are quite 
different from those expected for the interaction of ultra-high 
energy neutrinos. However, if a simulation starting from basic 
principles allows for a good reproduction of the experimental 
results, it is reasonable to assume that these simulation methods 
are transferable to neutrino interactions, as they are governed by 
the same underlying physical processes. 



beam 

(a) Top view 



(b) Isometric view 


3. Experimental Setup and Beam Characteristics 

The experiments presented in this paper were performed with 
a pulsed infrared Nd;YAG laser facility (A = 1064 nm) located 
at the Erlangen Centre for Astroparticle Physics (ECAP) of the 
University of Erlangen, and the pulsed 177 MeV proton beam 
of the “Gustaf Werner Cyclotron” at the “Theodor Svedberg 
Laboratory” in Uppsala, Sweden. The beam properties allow 
for a compact experimental setup. In both cases, the beams 
were dumped into a dedicated 150 x 60 x 60 cm^ water tank, 
where the acoustic held was measured with several position- 
adjustable acoustic sensors (see Eig.[T]). The sensors (also called 
hydrophones) could be positioned within the tank with absolute 
uncertainties below 1 cm. The temperature of the water could 
be varied between l°Cand 20°C with a precision of0.1°C. The 
temperature was brought to a particular value by hrst cooling 
the water with ice; subsequently the whole water volume was 
heated to the desired temperature in a controlled, gradual pro¬ 
cedure. Once the water temperature had been established, at 
least 10 min remained for measurements until the water volume 
heated up by 0.1°C through heat transfer from the environment. 
This time span was sufficient for all measurements conducted 
at water temperatures below the ambient temperature. 

The explored range of spill energies for the proton beam was 
from 10 PeV to 410 PeV, the beam diameter was approximately 
3 mm and the spill time U 30 p.s. Eor 177 MeV protons, 
the energy deposition in the water along the beam axis (z-axis. 


Figure 1: Schematic view of the experimental setup. Beams can enter the water 
volume through a thin foil (proton beam) or glass window (laser beam), which 
are set away from the tank walls via an air-filled tube to minimise effects of 
surface reflections on the recorded signals. Either of the two tubes, labelled ‘1’ 
and ‘2’ in |(a)| was used, depending on the experimental setup. The origin of 
the coordinate system was chosen at the entry point of the beam into the water 
target, see example in|(a)| 


beam entry into the water at z = 0 cm) is relatively uniform up 
to z = 20 cm ending in the prominent Bragg-peak at z ~ 22 cm 
(see Pig. 0. To adjust the spill energy, the number of protons 
per bunch was varied. The total charge of a bunch was cali¬ 
brated with two independent methods (Paraday cups and scintil¬ 
lation counters), leading to an uncertainty on the order of 15%, 
with some higher values for low spill energies. To obtain the 
beam intensity and prohle for the proton interactions in the wa¬ 
ter tank, the distance of about 1.2 m that the beam was travelling 
from the exit of the beam pipe through air and its entering into 
the water tank were included in the GEANT4 simulation. 

Eor the laser experiment, the pulse energy was adjusted be¬ 
tween 0.1 EeV and lOEeV and calibrated using a commer¬ 
cial power meter. The beam had a diameter of approximately 
2 mm and the pulse length was hxed at 9 ns. Eor the infrared 
light used, the laser energy density deposited along the beam 
axis has an exponential decrease with an absorption length of 
(5.9 ± 0.1) cm (see Pig.0. 

Eor both beam types the lateral energy deposition prohle was 
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Gaussian (the aforementioned beam diameters are the cr’s of the 
profiles). The two experiments allow to explore different spatial 
and temporal distributions of the energy deposition as well as 
two different mechanisms of energy transfer into the medium. 
For both beams, energy is deposited via excitation, in addition 
the medium is ionised in the case of the proton beam. 



Figure 2: Simulated energy deposition density along the beam axis (z-axis) for 
the two beam types (integral normalised to 1). For the proton beam the Bragg 
peak is prominent at 0.22 m, for the laser beam the decrease is exponential with 
an attenuation length of 0.059 m. 

For the signal recording, sensors based on the piezo-electric 
effect m were used. A full characterisation of these sensors 
had been performed prior to the experiments. They are linear 
in amplitude response, the frequency response is flat starting 
from a few kHz up to the main resonance at 40 kHz with a 
sensitivity of -155 dB re IV/qPa (~0.02 V/Pa). The main res¬ 
onance is more sensitive by ~5dBrelV/pPa and sensitivity 
drops rapidly at higher frequencies; at 90kHz, the sensitivity 
has dropped by 20dB. The absolute uncertainty in the determi¬ 
nation of the sensitivity is at a level of 2 dB in the frequency 
range of interest. Above 90 kHz the uncertainty exceeds 5dB. 
To calculate the response of the sensors to an external pres¬ 
sure pulse a parametrised fit of an equivalent circuit model as 
described in ifTSll was used. The sensitivity dependence on tem¬ 
perature was measured and the relative decrease was found to 
be less than 1.5% (or about 0.13 dB) per 1°C. 

For every set of fixed experimental parameters (temperature, 
energy, sensor position, etc.) the signals of 1000 beam pulses 
were recorded with a digital oscilloscope at a sampling rate in 
excess of 1 MHz. This rate is sufficient for the signals with 
spectral components up to 100 kHz, where the sensitivity of the 
sensors is negligibly small. These individual pulses were av¬ 
eraged to reduce background and environmental noise in the 
analysis, thereby obtaining a very high statistical precision. 

4. Basic Features of the Measured Signals 

Figurej^shows typical signals measured in the proton and the 
laser experiment, respectively, using the same sensors and ex¬ 
perimental setup. The general shapes of the two signals differ; 


a typical signal for the proton beam shows a bipolar signatur^ 
the one for the laser deviates from such a generic form. The 
laser signal has high frequency components up to several MHz 
due to the high energy deposition density at the point of beam 
entry and the almost instantaneous energy deposition compared 
to the ps pulse of the proton beam; therefore the resonance of 
the sensor is excited causing a ringing in the measured signal. 
The spatial distribution of the energy density e deposited by the 
laser leads to the two separate signals: the first originates in the 
beam area at the same z-region as the sensor placement (‘direct 
signal ’), the second from the beam entry, a point of discontinu¬ 
ity where most of the energy is deposed (‘beam entry signal’). 

The signal of the proton beam is deteriorated with respect to 
an ideal bipolar signal. Three main contributions to this dis¬ 
tortion can be discerned; the recorded signal starts before the 
expected onset of the acoustic signal (55.2 ps for the given po¬ 
sition, see Fig. given by the sonic path length); reflections 
of the acoustic wave on the beam enPy window overlay the 
original wave starting in the first rare-faction peak; and Anally 
there are frequency components of the signal exciting a reso¬ 
nant response of the sensor, slightly changing the signal shape 
and causing ringing. The first point was studied and found to be 
consistent with an elecPic charge effect in the sensors caused by 
the proton beam. Its starting time was always coincident with 
the beam pulse entry into the water, even for sensor distances 
of up to 1 m, hinting at an electromagnetic origin of the distort¬ 
ing signal. The shape of this non-acoustic signal is consistent 
with the integrated time-profile of the beam pulse with a subse¬ 
quent exponential decay. This deformation of the main signal 
is considered a systematic uncertainty on the signal properties 
and treated as such in the analysis. For the most part, its shape 
was fitted and subtracted from the signal. 

In order to minimise the impact of the signal deformation 
caused by the described effects on the analysis of the recorded 
signals, robust characteristics were used: the peak-to-peak am¬ 
plitude and the signal length from maximum to minimum of the 
signal. For the laser experiment these features were extracted 
for the direct signals only. 

5. Simulation of Thermo-Acoustic Signals 

For an in-depth validation of the thermo-acoustic model, 
comparisons of the signal properties with simulation results 
based on the model are essential. To this end, a simulation 
of the expected signals was developed. It is based on the 
thermo-acoustic model using a numeric solution of Eq. 

The input parameters to the simulation were either measured 
at the experiments, i.e. medium temperature and beam pro¬ 
files, or simulated, i.e. the energy deposition of the protons (us¬ 
ing GEANT4). The thermodynamic parameters bulk volume 
expansion coefficient, heat capacity and speed of sound were 
derived from the measured water temperature using standard 
parametrisations. Tap water quality was assumed. 


^Note that also the expected signal shape for acoustic pulses produced in 
interactions of ultra-high energy neutrinos in water is a bipolar one E). 
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Figure 3: Typical signals measured in the proton (upper graph) and laser (lower 
graph) experiment. Both signals (solid lines) were taken at a sensor position 
near .Y=10cm and z=20cm and are shown within the same time interval. The 
dash-dotted line in the upper graph indicates the charge effect described in the 
text. 


A series of simulations was conducted, where the input pa¬ 
rameters were varied individually in the range given by the ex¬ 
periment, including uncertainties. Especially the spatial and 
temporal beam profile have a substantial impact on amplitude, 
duration and shape of the signal. Simulated signals and the 
respective sensor response corresponding to the measured sig¬ 
nals of Fig. are shown in Fig. To minimise systematic ef¬ 
fects from the setup caused e.g. by reflections on the surfaces, 
the sensor response was convoluted onto the simulated signals, 
rather than deconvoluted from the measured ones. Thus in the 
analysis voltage rather than pressure signals are compared. 

The shapes of the simulated pulses are altered by the sen¬ 
sor response, especially the high-frequency components above 
the resonant frequency of the sensor. In the case of the laser 
pulse, mostly the resonance of the sensor is excited, leading to 
a strong ringing. For the proton beam, the primarily bipolar 
shape is again prominent, whereas the laser pulse is segmented 
into the two parts described above. The direct pulse of the laser 
experiment exhibits a bipolar shape as well, albeit less symmet¬ 
ric than for the proton beam. 

Figure exemplifies the dependency of the signal ampli¬ 
tude on the input parameters of the simulation; water temper¬ 
ature, pulse energy, beam profile in x and y, pulse length and 
the position of the sensor. As nominal positions of the sensor 
X = 0.10 m, y = 0.0 m, and z = 0.22 m were used. All param¬ 
eters were varied by +50% around the value of the best agree¬ 
ment with measurement, i.e. the values used for the simulations 
of the signals shown in Fig. 

Some of the characteristics of thermo-acoustic sound gener¬ 


Figure 4: Simulated signals for the proton (upper graph) and laser experiment 
(lower graph). Both signals (solid lines) were calculated for the experimental 
conditions of Fig. The corresponding dashed-lined signals mark the simu¬ 
lated signals after convolution with the sensor response. For better comparabil¬ 
ity, the signal maxima were normalised to 1. 


ation are observable. As discussed in Sec.|^ the dependence on 
temperature enters through the factor v^afcp (where the depen¬ 
dence on the speed of sound is negligible) and is roughly linear 
in the range investigated for this study. The dependence on en¬ 
ergy is strictly linear. The dependency on the beam pulse pa¬ 
rameters is diverse. It is governed by the integral in Eq. ( [T0| and 
therefore depends on both the spatial and the temporal beam 
profiles and the interaction of the particles with the water. For a 
given point in space and time, the elementary waves produced 
in the volume of energy deposition may interfere constructively 
or destructively depending on the beam properties. Accord¬ 
ingly, the length of the laser pulse has no influence on the am¬ 
plitude, as with 9 ns it is much shorter than the transit time of 
the acoustic signal through the energy deposition area. For the 
several ten ps long proton spill, the spill time is comparable to 
the transit time. Thus the acoustic signal shows a strong depen¬ 
dence on the spill time. The dependence on the radial coordi¬ 
nate r w.r.t. the beam axis r follows roughly the expected Ij ^|r 
and 1 /r fall-off of a cylindric source in the near and the far held, 
respectively. The y-position was varied between -2 and 2 cm, 
as the signals were recorded within the +z-plane. The resulting 
change in amplitude is below 1%. The z-dependency for the 
laser experiment follows the exponential fall-off expected from 
the light absorption. The one for the proton experiment shows 
the only non-strictly monotonic behaviour due to the form of 
the energy deposition with the prominent Bragg peak. 

Using these dependencies the systematic uncertainties of the 
model were obtained using the experimental uncertainties of 
the various parameters. The main uncertainty for the proton 
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rel. input variation 

(a) Proton beam simulation results 



rel. input variation 

(b) Laser beam simulation results 


Figure 5: Dependency of the simulated signal amplitude on variation of the 
simulation input parameters (labelled a - h) for the proton (upper graph) and 
the laser experiment (lower graph). With a nominal position y = Om of the 
sensor, the full range of the variation corresponds to ±2 cm for this parameter. 


beam is given by the temporal profile of the pulse, which was 
simplified to a Gaussian profile for this chapter (however not 
for the rest of this work). For the laser beam experiment this 
parameter influences the signal amplitude only on a one percent 
level. The second main influence on the amplitude is the sensor 
position along the beam axis (z-direction). 


Table [T] gives the parameters and their uncertainties (AV) 
used for the simulation of the signals. The resulting system¬ 
atic uncertainties in the amplitude (AA) are given as well. The 
combined uncertainties are ^ 29 % for the proton signal and ^^ 4 % 
for the laser signal, respectively. 


Parameter Proton Laser 



value (AV) 

AA 

value (AV) 

AA 

temperature 

12.0°C 

(+ 8 %) 

+ 10% 
-11% 

12.0°C 

(+ 8 %) 

+ 11 % 

energy 

410PeV 

(+ 10 %) 

+ 10 % 

520 PeV 
(+ 10 %) 

+ 10 % 

beam width 
(x-direction) 

3.0 mm 
(+14%) 

-14% 

+ 19% 

5.0 mm 
(+15%) 

+ 1 

beam width 
(y-direction) 

3.0 mm 
(+14%) 

-12% 

+ 16% 

5.0 mm 
(+15%) 

+ 0 % 

pulse length 

1 . 0 - 10 -^ s 
(+ 11 %) 

-13% 

+ 16% 

6.9- lO-’s 
(+50%) 

-2% 

+2% 

x-position 

y-position 

0 . 10 m 

(+7%) 

0.0 m 
(+1 cm) 

-4% 

+5% 

-2% 

+0% 

0 . 11 m 

(+ 6 %) 

0.0 m 
(+1 cm) 

+4% 

-1% 

-1% 

z-position 

0.22 m 
(+4%) 

+9% 

0.17m 

(+5%) 

-11% 

+ 13% 

combined 


+34% 

-29% 


+26% 

-24% 

Table 1: Beam parameters used for the simulated signals in Fig. with their 
associated experimental uncertainties (AV) and resulting uncertainty in signal 
amplitude (AA). The pulse length of the laser is set to a value much higher 


than in the experiment, to save calculation time. The resulting uncertainty in 
the signal amplitude is below 2% 


6. Analysis Results 

In the following, the results obtained in the comparison of 
experiment and simulation are presented. 

6.1. Comparison of Simulated and Measured Signals 

6.1.1. Proton Beam Experiment 

Figure shows a comparison between simulated and mea¬ 
sured signals for the proton beam experiment at different sensor 
positions. 

For better visibility only the main part of the signals (first 
bipolar part) is plotted. The input parameters of the simulation 
were varied within the experimental uncertainties until the best 
agreement with the measured signal in amplitude and duration 
was obtained for the reference point at x=0.40 m and z-Q. 11m. 
For this optimisation, a simple procedure of adjusting the pa¬ 
rameters manually and scanning the resulting agreement visu¬ 
ally was found to be sufficient. The other signals were sim¬ 
ulated with the same parameter set, only the sensor positions 
were changed. The signal shapes differ for different z-positions 
due to the geometry of the energy deposition profile described 
in Sec. |^with cylindric form in the xy-plane and almost flat en¬ 
ergy density in z-direction up to the Bragg peak at z-Q.22 m. 
Due to this geometry an almost cylindrical wave is excited in 
the medium, with coherent emission perpendicular to the beam 
axis. In the region z < 0.2 m the signals are of bipolar shape. 
Along the beam axis (x 0.0 m, z > 0.25 m) the main part of 
the observed signal originates from the Bragg peak as a nearly 
spheric source and no clear bipolar shape evolves. 
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(c) Signal at ;i: = 0.01 m, z = 0.41 m 


Figure 6: Comparison of measured (solid line) and simulated (dashed line) 
signals at different sensor positions for the proton beam experiment. Note the 
different scale of the x-axis used in plot |(c)| The reference signal for tuning the 
simulations is at z = 0.11 m and x = 0.40 m. The coordinate system is given in 
Fig.^ See text for a detailed description. 



Figure 7: Comparison of measured (solid line) and simulated (dashed line) sig¬ 
nals for the laser beam experiment at z = 0.17 m for different sensor positions 
along X (reference at x = 0.10 m). For more details see text. 


The agreement between simulation and measurement is good 
for all positions. Not only amplitude and duration match (see 
also the following sections) but also the signal shape is repro¬ 
duced to a very high degree. The small discrepancies, primarily 
in the rare-faction part of the bipolar pulse, have contributions 
stemming from a non-ideal sensor calibration and reflections on 
the tank surfaces. However, a significant part of the discrepan¬ 
cies may lie in the beam pulse modelling or even the thermo¬ 
acoustic model itself. 

6.1.2. Laser Beam Experiment 

The prominent feature of the energy deposition of the laser 
beam is at the beam entry into the medium. Overall, the ge¬ 
ometry of this deposition is mainly a cylindric one with rota¬ 
tional symmetry around the z-axis, as for the proton beam, with 
coherent emission perpendicular to the beam axis (direct sig¬ 
nal). The signal from the discontinuity at the beam entry is 
emitted almost as from a point source (beam entry signal). In 
contrast to the proton beam, the shapes of signals at different 
positions along the z-axis do not vary much, only the relative 
timing between the two signal components varies. Therefore, 
Fig. 1^ shows only signals for a z-position in the middle of the 
water tank. Here, both signal parts are well described by the 
simulation. The beam entry part of each signal is less well re¬ 
produced in the simulation due to its high frequency compo¬ 
nents where the sensor calibration is less well understood. 

To compare the characteristics of the signals in simulation 
and measurement in the following studies, their amplitudes and 
the point in time of the recording of the signal maximum were 
further studied. For the laser signal, the direct signal part is 
considered only. 

6.2. Speed of Sound Measurement 

Equation yields as velocity of propagation for the 
thermo-acoustic signal the speed of sound in the medium, here 
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water. To verify the hydrodynamic origin of the measured sig¬ 
nals the variation of the arrival times for different sensor po¬ 
sitions perpendicular to the beam axis were analysed. Fig¬ 
ure in shows the measured data and a linear fit for each beam 
type. The data is compatible with an acoustic sound propa¬ 
gation in water, as the fits yield a speed of sound compatible 
with pure water at the temperature used. For the proton beam 
Vj = 1458 + 4 m/s (j^f^/n.d.f. = 8.2/11) was obtained for a 
water temperature of 11.0 + 0.2 °C, where literature M gives 
for pure water at normal pressure = 1455 + 1 m/s which is 
in complete agreemenj^ For the laser measurements the wa¬ 
ter temperature of 19.4 + 1.0 °C and thus the speed of sound 
were significantly higher, the observed Vj = 1508 + 3 m/s 
(/p^/n.d.f. = 9.2/13) is again in perfect agreement with the the¬ 
oretical value of 1503 + 6 m/s. The offset between proton and 
laser beam data in Fig. is due to a differing delay time in be¬ 
tween trigger time and arrival time of the different beams in the 
water and is irrelevant for the calculation of the speed of sound. 



Figure 8: Arrival time of the signal maximum for different sensor position per¬ 
pendicular to the beam (x position). The straight lines represent linear fits to 
the data points yielding v, = 1458 ± 4 m/s and v,, = 1508 ± 3 m/s for the proton 
and laser beam, respectively. 


6.3. Variation of the Sensor Position 

The durations of the complete, undipped signals vary with 
sensor positions mainly due to reflections, which were not sim¬ 
ulated and can therefore not be compared. 

The comparison of simulated and measured signals for differ¬ 
ent sensor positions within the water tank, excluding the parts of 
the signals dominated by reflections, was shown in Fig. The 
good agreement between model and measurement also mani¬ 
fests itself in the development of the signal amplitude with dis¬ 
tance of sensors from the beam axis, shown in Fig. To min¬ 
imise systematic effects from reflections, only the amplitude of 
the leading maximum is analysed. 


^The uncertainties of the theoretical values of the speed of sound result from 
the changing water temperature while the measurements were taken. 


Though there are sizeable deviations, the overall shape of the 
curve is reproduced. The behaviour is again different for the 
two z-positions. The development at z = 0.11m follows the 
one expected from a cylindric source with alj xjr behaviour in 
the near field up to ~ 0.3 m and a 1/r behaviour in the far field 
beyond that distance. At the smallest measured distances, the 
simulated behaviour deviates from the measured one. Presum¬ 
ably this is due to simplifications made in the derivation of the 
model in Sec.|^ At z = 0.21 m the amplitude falls off more uni¬ 
formly, this is again a combination of the point-like emission 
characteristic of the Bragg peak interfering with the cylindric 
emission at z < 0.2 m. 



Figure 9: Development of the leading peak (maximum) with sensor position 
perpendicular to the beam axis (x) for the proton beam experiment. Points 
mark the measured amplitudes for two different z-positions and the lines the 
respective simulation results. 


The behaviour for the laser experiment is not as well repro¬ 
duced by the simulation (see Fig. 10 1 . This is again attributed to 
the high-frequent signal components of the laser, where minor 
uncertainties in the simulation may result in big changes of the 
signal amplitude. Especially the overlap of the direct with the 
beam entry window signal distort the signal shape. At distances 
exceeding 0.5 m the two signal parts cannot be distinguished. 


6.4. Variation of the Pulse Energy 


Assuming otherwise unchanged settings, the energy deposi¬ 
tion density e scales linearly with total deposited energy. Thus 
the spill energy can be written as a pre-factor in Eq. ( [TOj i effect¬ 
ing the pressure and thus signal amplitude linearly. As shown in 


Eig. 11 this behaviour was observed in the experiments yielding 
a zero-crossing of the pulse energy at (-4.1 + 5.3) mPa for the 
proton beam and (42 ± 87) mPa for the laser beam. Both values 
are consistent with zero. The slope of the line depends on the 
energy deposition and the sensor positioning along the beam 
axis and can therefore not be compared for the two beams. As 
expected from the model, the signal duration and signal shape 
showed no significant dependence on energy. 
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Figure 10: Development of the leading peak (maximum) with sensor position 
perpendicular to the beam axis (v) for the laser beam experiment. Points mark 
the measured amplitudes at z = 0.27 m and the line the simulation result. 


6.5. Variation of the Temperature 


The main feature of the thermo-acoustic model is its depen¬ 
dence on the temperature of the medium. Figure shows the 
temperature dependence of the signal peak-to-peak amplitude 
for the laser beam, where a positive (negative) sign denotes a 
leading positive (negative) peak of the signal. The two data sets 
shown in the figure were recorded by two sensors simultane¬ 
ously, which were positioned at x = 0.10 m perpendicular to 
the beam axis and at z = 0.11 m and z = 0.21 m along the beam 
axis, respectively. In the case of the proton beam setup, which 
will be discussed below, these hydrophone positions correspond 
roughly to the z-position of the Bragg-peak and a z-position half 
way between the Bragg-peak and the beam entry into the wa¬ 
ter, respectively. For comparability, the same positions and the 
same sensors were chosen for the laser and proton beam exper¬ 
iments. 

The laser beam signal shown in Fig. [T^ changes its polarity 
around 4°C, as expected from the thermo-acoustic model. The 
theoretical expectation for the signal amplitude, which is pro¬ 
portional to ajcp and vanishes at 4°C for the given temperature 
and pressure, is fitted to the experimental data. In the fit, an 
overall scaling factor and a shift in temperature (for the exper¬ 
imental uncertainty in the temperature measurement) were left 
free as fit parameters. The fit yielded a zero-crossing of the am¬ 
plitude at (3.9 + 0.1)°C, where the error is dominated by the 
systematic uncertainty in the temperature setting. 

Analysing the proton data in the same fashion resulted in a fit 
that deviated from the model expectation, and a zero-crossing 
significantly different from 4.0°C at (4.5 + 0.1)°C, see Fig. 13 
The data strongly indicate the presence of a systematic effect 
near the zero-crossing of the signal amplitude. To understand 
this effect, the signal shapes near the temperature of 4.0°C were 
investigated (Fig.[T4|). A non-vanishing signal is clearly observ¬ 
able at 4.0°C and the signal inverts its polarity between 4.0°C 
and 4.5°C. In view of the results from the laser beam measure¬ 
ments and the obvious systematic nature of the deviation from 



Figure 11: Peak-to-peak amplitude of signals for dilferent pulse energies, for 
both beams the sensors where at a position x « 0.10 m. There is a ten percent 
systematic uncertainty in the absolute determination of the pulse energy. The 
lines represent linear fits to the data points yielding a zero-crossing of the am¬ 
plitude compatible with no energy in a pulse. The insert shows the data for the 
proton beam. 


the model visible in Fig. we subtracted the residual signal 
at 4.0°C, which has an amplitude of 5% of the 15.0°C signal, 
from all signals. Thus a non-temperature dependent effect in 
addition to the thermo-acoustic signal was assumed. The re¬ 
sulting amplitudes shown in Fig. 15 are well described by the 
model prediction. 


The production mechanism of the underlying signal at 4.0°C, 
which was only observed in the proton experiment, could not be 
unambiguously determined with the performed measurements. 
From the model point of view, the main simplification for the 
derivation of Eq. ( [T0| was to neglect all non-isotropic terms 
and momentum transfer to the medium in the momentum den¬ 
sity tensor IT? by setting /? = 0 in Eq. (j^. As discussed in 
Sec. dipole radiation could contribute significantly near the 
disappearance of the volume expansion coefficient for the case 
P + 0. Also other non-thermo-acoustic signal production mech¬ 
anism have been discussed in the literature which could give 
rise to an almost temperature independent signal, see e.g. Q. 
The obvious difference to the laser experiment are the charges 
involved both from the protons themselves and the ionisation 
of the water which could lead to an interaction with the polar 
water molecules. Another difference are the massive protons 
compared to massless photons. Residual signals at 4.0°C were 
found in previous works as well 1^61 [H, as will be discussed 
in more detail in Sec. 16.61 


Eor clarification further experiments are needed either with 
ionising neutral particles (e.g. synchrotron radiation) or with 
charged particles (e.g. protons, a-particles) with more sensors 
positioned around the Bragg-peak. With such experiments it 
might be possible to distinguish between the effect of ionisation 
in the water and of net charge introduced by charged particles. 


9 


























































Figure 12: Measured signal amplitude of the bipolar acoustic signal induced by 
laser pulses at different temperatures fitted with the theoretical expectation as 
described in the text. All amplitudes were normalised to 1 at 15.0°C. The insert 
shows a blow-up of the region around 4° C where the sign of the amplitude 
changes. 


6 .6. Comparison with Previous Experiments 

With the analysis that has been described above, the sig¬ 
nal production according to the thermo-acoustic model could 
be unambiguously confirmed. While the model has been con¬ 
firmed in previous experiments, the simulations presented in 
this work constitute a new level of precision. The most puz¬ 
zling feature, a residual signal at 4°C that was also observed 
in previous experiments, was investigated with high precision 
by scanning the relevant temperature region in steps of 0.1 °C. 
The observed shift of the zero-crossing of the amplitude to¬ 
wards values higher than 4°C, caused by a leading rarefaction 
non-thermal residual signal at 4°C, is in qualitative agreement 
with Q- In IIH, a residual signal at4°C was also reported. Since 
in that work the zero-crossing of the amplitude is observed at 
6 °C, i.e. a higher value than the expected 4°C, it can be assumed 
that the corresponding residual signal has a leading rarefaction. 

In 0, a residual signal is found at 4.25°C, however with a 
leading compression rather than rarefaction. The authors con¬ 
clude that this may lead to a signal disappearance point below 
the expected value, in contrast to |l4l|5l and the work presented 
in this article. 

For the measurements with a laser beam reported in Q, a 
residual signal was also observed at 4°C, albeit with a leading 
compression and a subsequent reduction of the temperature of 
the zero-crossing of the signal amplitude to about 2.5° - 3.0°C. 
This observation is in contrast to the laser experiment presented 
in this article. 

In conclusion, the works of all authors discussed here indi¬ 
cate a non-thermal residual signal for proton beams, albeit with 
varying results concerning the size of the effect and the shape 
of the underlying non-thermal signal. The results for the laser 
beam reported in differ from those described in the article 
at hands. It should be pointed out, however, that in ||5l and 
0 results are reported by the same authors for proton and for 



Figure 13: Measured signal amplitude of the bipolar acoustic signal induced 
by proton pulses at dilferent temperatures near 4.0° C, fitted with the model 
expectation as described in the text. A systematic deviation from the model 
expectation as the amplitude changes its sign is clearly visible. The amplitudes 
were normalised to 1 at 15.0°C. 


laser beams, respectively. A comparison of these two publica¬ 
tions shows that a different behaviour near the temperature of 
4°C was observed for the two types of beams. Hence, to the 
best of our knowledge, there are currently no results in contra¬ 
diction with the notion of different non-thermal effects in the 
interaction of proton vs. laser beams with water. The available 
data does not allow for a more detailed analysis of the corre¬ 
lation between experimental conditions and the temperature of 
the zero-crossing of the signal amplitude. 

7. Applications in Astroparticle Physics 

Efforts to detect neutrinos at ultra-high energies are at the 
frontier of research in the field of astroparticle physics. Neu¬ 
trinos are the only viable messengers at ultra-high energies be¬ 
yond the local universe, i.e. distances well beyond several tens 
of Megaparsecs. If successful, the investigation of these elusive 
particles will not only enhance the understanding of their own 
nature, but also provide important complementary information 
on the astrophysical phenomena and the environments that ac¬ 
celerate particles to such extreme energies. 

For acoustic particle detection, not only the technical aspects 
such as optimal design and detector layout are subject of re¬ 
search. But also the underlying physics processes—i.e. the 
formation of hadronic cascades resulting from neutrino inter¬ 
actions in dense media—have never been observed directly in 
detector experiments at these energies. Producing ever more re¬ 
liable extrapolation of reaction properties to ultra-high energies 
is an ongoing effort. With advancements in the simulations of 
cascades forming in water and improvements of detector simu¬ 
lation tools, the discrepancy between cascade parameters from 
independent simulations decreased: recent studies differ only 
slightly iniiisMi. At the same time it is necessary to gain 
a solid understanding of the sound signals generated from the 
energy depositions by particle cascades. For this purpose, lab¬ 
oratory measurements are required. 
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Figure 14: Comparison of signals recorded for the proton beam with hy¬ 
drophone 1 at different temperatures. For each temperature, the signals were 
first smoothed with a 100 kHz second order Butterworth low pass filter. Subse¬ 
quently 1000 waveforms were averaged, as described in Sec.[^ To allow for an 
easy comparison of the signal shapes, a point in time at the onset of the acoustic 
signals was chosen as zero time and for all signal amplitudes the corresponding 
offset was added or subtracted to yield a zero amplitude at that time. 


This work, together with others ll4U8l l^l2^ has established 
the validity of the thermo-acoustic model with uncertainties at 
the 10% level. The input to a model as discussed in this article is 
an energy deposition in water as it is also produced by a cascade 
that evolves from a neutrino interaction. In comparison with 
the uncertainty in the thermo-acoustic model, uncertainties due 
to the simulations of hadronic cascades and cascade-to-cascade 
variations are large im, dominating the challenge to detect and 
identify sound signals resulting from neutrino interactions. It 
can hence be concluded that the current level of precision in 
modelling sound signals in the context of the thermo-acoustic 
model is fully sufficient for the understanding of acoustic neu¬ 
trino signatures. The latter is necessary to improve the selec¬ 
tion efficiency and background rejection for neutrino detection 
algorithms in potential future acoustic neutrino detectors. Sev¬ 
eral experiments have been conducted lISTl 12414^ to under¬ 
stand the acoustic background at the sites of potential future 
large-scale acoustic neutrino telescopes in sea water, fresh wa¬ 
ter and ice. The combination of simulation efforts, laboratory 
measurements and studies with in-situ test arrays will allow for 
a conclusion of the feasibility of acoustic neutrino detection. 

8. Conclusions 

We have demonstrated that the sound generation mecha¬ 
nism of intense pulsed beams is well described by the thermo¬ 
acoustic model. In almost all aspects investigated, the signal 
properties are consistent with the model. The biggest uncer¬ 
tainties of the experiments are on the 10% level. One discrep¬ 
ancy is the non-vanishing signal at 4°C for the proton beam 
experiment, which can be described with an additional non¬ 
temperature dependent signal with a 5% contribution to the am¬ 
plitude at 15°C. The model allows for calculations of the char¬ 
acteristics of sound pulses generated in the interaction of high 



Figure 15: Measured signal amplitude of the bipolar acoustic signal induced 
by proton pulses at different temperatures fitted with the model expectation as 
described in the text. The non-thermo-acoustic signal at 4.0° C was subtracted at 
every temperature. The amplitudes were afterwards normalised to 1 at 15.0°C. 
The insert shows a blow-up of the region around 4° C where the sign of the 
amplitude changes. 


energy particles in water with the input of the energy deposi¬ 
tion of the resulting cascade. A possible application of this 
technique would be the detection of neutrinos with energies 
> 1 EeV. 
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